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Abstract 

We extend the SafFman theory of membrane hydrodynamics to account for the correlated motion 
of membrane proteins, along with the effect of protein concentration on that correlation and on 
the response of the membrane to stresses. Expressions for the coupling diffusion coefficients of 
protein pairs and their concentration dependence are derived in the limit of small protein size 
relative to the inter-protein separation. The additional role of membrane viscosity as determining 
the characteristic length scale for membrane response leads to unusual concentration effects at large 
separation — the transverse coupling increases with protein concentration, whereas the longitudinal 
one becomes concentration- independent. 

Key words: membrane hydrodynamics; Brownian motion; hydrodynamic interaction; effective vis- 
cosity; diffusion coefficient 



Introduction 

Biomembranes contain a high concentration of proteins, performing key cellular functions (jH). Ex- 
tensive efforts have been directed, therefore, at measuring the dynamics of membrane proteins us- 
ing various experimental techniques (Q)- Those studies have concentrated on either single-molecule 
dynamics (using single-molecule tracking (0, 0) and fluorescence correlation spectroscopy (j2])) or 
collective gradient diffusion (using fluorescence recovery after photobleaching (jsl)). Both the single- 
particle and large-scale levels are in some contrast with the expected cooperative motion of several 
proteins at the high concentration relevant to biomembranes, and the crucial role played by small 
protein aggregates in membranes. Until now experimental investigations of the correlated motion 
in membranes and monolayers have been limited to rather extended objects, such as domains ([lj0) 
and embedded microspheres (@), although similar two-point microrheological measurements with 
fluorescent membrane proteins seem feasible (0)- Furthermore, in cases where one is interested in 
membrane properties rather than those of the proteins, two-point microrheology has the advantage 
of being insensitive to the shape of the inclusion and the local perturbation that it introduces in 
the membrane (|9l. IToh. It seems to be of significant interest, therefore, to account for the correlated 
Brownian motion of membrane proteins and the associated effects on membrane dynamics. 

From a hydrodynamic perspective, a bare, protein-free membrane can be viewed as a quasi-two- 
dimensional (quasi-2D) liquid, whose molecules (lipids) are free to flow only within the membrane 
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surface, yet exchange momentum with the surrounding three-dimensional (3D) hquid. Hence, flows 
within a membrane are essentially different from those in both 3D and 2D liquids, in that they 
do not conserve momentum. As a result, a characteristic length scale emerges, such that over 
distances much smaller than the membrane keeps its momentum and responds similar to a 
2D liquid, whereas beyond that distance momentum is exchanged with the surroundings, and the 
response is significantly modified. 

To gain further intuition for the results that will follow, let us compare this situation in slightly 
more detail with the ones in momentum-conserving 3D and 2D liquids. In the 3D case the stress 
(i.e., momentum flux) emanating from a local perturbation must decay with distance r as u ^ 
(so that the total flux through an envelope of radius r should remain constant). Since shear stress 
is related to fluid velocity v a.s a ^ ijfVv, where rji is the shear viscosity of the fluid, the velocity 
response to the perturbation decays as w '--^ {r]{r)~^. More specifically, given a point force f acting 
on the liquid at the origin, the steady-state flow velocity v(r) at position r is given by (jlir i 



3D liquid: Wi(r) = (r)/j. 



G.,(r) 



87r77fr 



(1) 



where G is the Oseen tensor and i,j = x,y,z. (Summation over repeated indices is implied 
throughout this article.) An important consequence of this velocity response is that the correlation 
(hydrodynamic interaction) between the motions of two particles suspended in the liquid is long- 
ranged, decaying as 1/r ([ill). Similarly, in a 2D momentum-conserving liquid the stress must decay 
as (T ~ r~^, and, thus, the velocity decays as —f]^^ In(K'r), where ?7m is the 2D shear viscosity. The 
resulting well known logarithmic divergence of this problem (|ll[ ) requires a cutoff length k'~^ (e.g., 
the system size). The analogue of Eq. [T]for the 2D case is 



2D liquid: ^;,(r) = G„ (r)/,, G„ (r) 



1 



-ln(«;'r)% + l^ 



(2) 



where here (and in the rest of the article) i,j — x,y. For the intermediate quasi-2D case of fluid 
membranes we expect the response, and thus the hydrodynamic interaction between inclusions, to 
cross over from the 2D logarithmic decay at short distances to the 3D 1/r decay at large distances. 

Three major theoretical approaches to membrane hydrodynamics have been presented. In 
Saffman's pioneering theory (|l2l ) the membrane is modeled as a flat slab of viscous liquid, having 
width w and viscosity rjai/w. The flow velocity in the membrane is assumed to be two-dimensional, 
i.e., the velocity proflle across the slab width is uniform. The slab has at its two bounding sur- 
faces no-slip contacts with an inflnite viscous fluid (water), having viscosity rjf. The emergent 
characteristic length — the Saffman-Delbriick length (|l3l ) — is given by 



2r/f' 



(3) 



(In asymmetric cases, where the liquids on the two sides of the membrane have different viscosities, 
one should replace 27yf with the sum of those viscosities.) Considering a membrane protein as a 
cylindrical inclusion of radius a k~^, the Saffman theory yields the following self- mobility for 
the protein, 

B, = -^[-ln(W2)-7], (4) 

where 7 ~ 0.58 is Euler's constant. The self-diffusion coefficient of the protein is simply given, 
through Einstein's relation, by Dg — k^TB^, where k^T is the thermal energy. The weak logarith- 
mic dependence of Ds on the protein size, as predicted by Eq. |4l does not seem to be obeyed in 
practice |l3 ) , possibly because of the local disturbance that the protein creates in the membrane 
(jlOl). Equation |4] has been confirmed, nonetheless, in the Brownian motion of small membrane 
domains The Saffman theory has been extended to arbitrary values of na (flsh and to the case 
of membranes supported on a liquid layer of finite thickness (|l6l). 
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In the second approach, introduced by Levine et al. (17, 3, the membrane is treated as a 
viscoelastic fihTi of vanishing thickness within an infinite viscous hquid, taking into account both 
in-plane and out-of-plane dynamics. The in-plane response of this model, in the limit of a purely 
viscous, incompressible film, coincides with that of the Saffman theory. 

The third theoretical approach considers the membrane as an effective 2D Brinkman liquid 
((20h . i.e., a liquid with an additional phenomenological term that makes it lose momentum over 



distances larger than a certain given value, k~ (|21l l22l |23| . \24 . 125f ). The mobility of a disk of 
radius a, as calculated from this theory, coincides with Eq. |3]in the limit rea <C 1. However, this 
approach is essentially different from the first two in that it breaks the translational symmetry 
along the membrane surface. Thus, while the theories of Saffman and Levine et al. conserve total 
momentum in 3D, allowing the surrounding liquid to impart momentum back to the membrane 
at large distances, in the Brinkman-like theory the momentum, once leaving the membrane, is 
lost. As a result, the large-distance response of this model is qualitatively different, decaying as 
(as required by mass conservation in 2D ((26[ )1 rather than 1/r (as resulting from momentum 
conservation in 3D). This approach, therefore, is appropriate at sufficiently short distances, or when 
translational symmetry is indeed broken, as in membranes supported on a solid substrate (3, 21). 



Since the lateral size of membrane proteins is typically at least one order of magnitude larger 
than that of the lipids, they can be considered as suspended in a continuous quasi-2D liquid, thus 
making the membrane a quasi-2D suspension. In analogy with ordinary suspensions, we expect 
that the average effect of many such mobile proteins will lead to a modified effective response of 
the membrane. For example, the presence of hard spheres in 3D suspensions leads to a modified 



effective viscosity, which is given, up to linear order in the volume fraction of spheres cf), by (|27f ) 



3D suspension: rj^^ = rjf ^1 -t- . (5) 

The analogous result for a 2D suspension of hard disks is ((28l ) 

2D suspension: r^^^f = r;,„ (1 + 2(j)) , (6) 

where here (j) is the area fraction of disks. As described above, membranes represent a more 
complicated intermediate between 2D and 3D liquids, and we expect, therefore, a more subtle, 
distance-dependent effect of inclusions on its response. 

In the current work we extend the Saffman theory of membrane hydrodynamics (fT3 ) to account 
for the correlated Brownian motion of protein pairs and the effect of protein concentration on 
that motion. The analysis is restricted to the limit of small protein size, kg <C 1. Since the 
viscosity of a lipid bilayer is typically 10"^-fold that of water (|29r ). is typically three orders of 
magnitude larger than the membrane thickness, i.e., of micron scale. Hence, the limit ko <^ 1 
should hold for any membrane protein, as originally assumed by Saffman and Delbriick. As regards 
the inter-protein distance, two regimes are addressed — intermediate distances, a <C r <C k^^, and 
the asymptotically far region, r 3> k~^. We avoid the region r ~ a, in which specific effects of 
protein shape and membrane distortion are expected to be important. Corrections to leading order 
in a/r are nonetheless derived. 

We begin with a presentation of several basic properties of Saffman's hydrodynamic theory, 
which are useful for our calculations. The resulting coupling diffusion coefficients of an isolated 
protein pair are subsequently described. We then proceed to derive the leading corrections to the 
membrane response, as well as the coupling diffusion coefficients, due to a low concentration of 
membrane proteins. Finally, we discuss the physical meaning of the results and their practical 
limitations. While writing this article we have learned of an independent study by Henle and 
Levine of the effective viscosity of membranes with mobile inclusions Isoj). The relation between 
their results and ours is addressed in the Discussion. 
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Theory 

Membrane hydrodynamics 

Consider a flat membrane, lying on the xy plane. First, we address the steady-state flow velocity of 
the membrane at position r, v(r), in response to a point force, f , applied at the origin and directed 
along the membrane plane. In Fourier space, v(q) = / dre^''' ''v(r), Saffman's analysis yields f|12il 

5.(q) = G,,(q)/,, G,,(q) = J—- (5,,-^), (7) 

where i,j — x,y. The tensor G is the membrane-analogue of the Oseen tensor of Eq. [TJ Inverting 
to real space, we get 



Kr 2 T:\Kry 

'iTj 



7r(Kr)^ 



(8) 



where Yn and are, respectively, Bessel functions of the second kind and Struve functions. At 
short distances, r <^ k"^, G reduces to 

r<K-i: G.,(r)~^-^{-[ln(Kr/2)+7 + l/2]% + ^} + 0(Kr). (9) 

This result coincides with the one for a 2D liquid, Eq. [2l with an appropriate definition of n' ~ k. 
At large distances, r ^ k^^, G tends to 

r»K-i: G,,— -^^ + 0^-2, (10) 

which shows the typical 1/r decay of 3D flows. Moreover, since rj^K = 2ri{, the large-distance 
response depends solely on the outer liquid viscosity and is independent of membrane viscosity. (In 
fact, up to a numerical prefactor, Eq. [10] could be readily derived by requiring that G decay as 1/r 
and obey 2D incompressibility, diGij =0.) 

Now suppose that, rather than being a point force, the force is applied to the membrane by 
a disk-like particle of finite radius a. To leading order in a/r, the membrane flow velocity is 
given by Vi{r) ~ Gij{r)fj. We are interested in the finite-size correction to this flow while still 
neglecting terms of order Ka (as we do throughout this work). The domain of interest, therefore, 
is a < r <C K^^. In this region the membrane behaves as a 2D liquid, following Eq. [51 Thus, 
the calculation reduces to a 2D Stokes problem of finding the flow away from a rigid disk, whose 
solution is 

a<r««-i: v,{r) ^ G[f {r) f„ gJ) (r) = G., (r) + (^5., - ^) , (11) 

where G is given by Eq. [5] 

Next we consider a membrane with a preexisting flow velocity v(r), and embed in it a circular 
disk of radius a moving with linear velocity U and angular velocity ft. We would like to find the 
force F, torque L, and force dipole (stresslet) S, which the inclusion exerts on the fluid membrane. 
For a sphere in an unbounded liquid the linear relations between (F, L, S), on the one hand, and 
(v, U, rj), on the other, are given by the first and second Faxen laws ([11). In Appendix A we derive 
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the (approximate) membrane-analogues of these laws, which are as follows: 



F 
L 



ln(Ka/2) +7 



1 



v+ -a^V^v-U 

4 



1 



Si/ 



1 + -a^V^ ) (V X v) - 20 
1 



(12) 
(13) 
(14) 



Applying the first relation, Eg. 1121 to a disk moving in an otherwise stationary membrane (v = 0), 
we recover the Saffman-Delbriick mobility (Eq. [3]). These relations are exact for a 2D liquid (where, 
in addition, the term proportional to V^(V x v) in Eg . [T51 vanishes^ . whereas for membranes their 
validity is restricted to sufhciently small particles, ^ 1. (See Appendix A for details.) 



Correlated diffusion 



Consider a pair of membrane proteins undergoing Brownian motion while being separated by the 
2D vector r. We consider a time period t, which is sufficiently short such that r can be assumed 
constant, yet sufhciently long to yield Brownian displacements linear in t. We further assume that 
r is much larger than the protein sizes (radii), ai and 02. The displacements of the two proteins 
during time t obey the following relations: 



lArfAr'^) 



where Arf is the displacement of particle a (a = 1, 2) along the axis 



(15) 

x,y). The diffusion 



tensor D^^^ characterizes both the self-diffusion of the particles {a = (3) and the coupling between 
them (a. ^ (3). We define the x axis, without loss of generality, along the line connecting the pair, 
i.e., r — rx. This choice leads, by symmetry, to D]^y = 0. The coupled diffusion is then fully 
characterized by two coefficients: a longitudinal coupling diffusion coefficient, -D£(r) = Z?^^(rx), 
and a transverse one, D'^ (r) = (rx) . The former is associated with the coupled Brownian motion 
of the pair along their connecting line, and the latter with the coupled motion perpendicular to 
that line. 

For the overdamped dynamics considered here the diffusion tensor in Eq.[T5]is simply related to 
a pair-mobility tensor via the Einstein relation, D"^^ = k^TB'^^^ . The mobility tensor Bl^{r) gives 
the change in velocity Vi of particle 2 located at r due to a unit force in the j direction applied to 
particle 1 at the origin. In the limit a/r —^ {a being the larger of ai, 02) this, in turn, is just the 
membrane-analogue of the Oseen tensor, Eq. [51 Hence, we readily identify. 



r»(a2«-i)i/3 . 
Dlir) ~ kBTG^rx) 

D'^{r) ^ kBTGyyirx) 



kBT 
kBT ' 



Hi{Kr) -Yi{Kr) 



4?7n 



Hoinr) 



Hiinr) 



i(yo(..)-r2(.0) + ^ 



(16) 



(The domain of validity stated here will be clarified below, when we address the effect of finite 
particle size.) This result coincides with the in-plane response functions derived by Levine and 
MacKintosh in the limit of a purely viscous, incompressible membrane |l7[). 



At very large inter-particle distances, r ^ k ^, Eg . [TBI reduces to 



Dlir) 



kBT 



kBT 



mr) 



kBT 



27r?7nj(Kr)-^ Sirrif'^r'^ 



(17) 



The longitudinal coupling between the two proteins decays asymptotically as 1/r and is independent 
of membrane viscosity. (It is identical, in fact, to the analogous coefficient in a 3D liquid.) The 
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transverse coupling decays only as l/r^, and, curiously, increases with membrane viscosity. When 
the inter-particle distance is smaller than yet still sufficiently large. Eg. If 61 becomes 



1/3 



fceT 



-ln(Kr/2) - 7± f/2 + (f =F f/3)Kr] , 



(18) 



where the upper (lower) signs correspond to the longitudinal (transverse) coefficient. 

Let us now examine the leading effect of finite particle sizes, oi and 02, and see at what inter- 
particle distance this effect becomes significant. In this domain, clearly, r ^ k~^. First, according 
to Eq. [TT] there is a correction of order af/r'^ to the flow velocity caused by the forced particle 
I. In addition, the first Faxen law, Eq. [121 yields a correction of order a^/r'^ for the velocity of 
particle 2 as it is embedded in that flow. Substituting v(r) of Eq. [TT]in Eg. [T^ while setting F = 
(particle 2 being force- free), we find the velocity U of particle 2 to leading order in a\,a\. The 
relation between U and f defines a corrected pair mobility tensor, resulting in the following coupling 
diffusion coefficients: 



-ln(Kr/2)-7±l/2±^ 



2r2 



(19) 



where the plus (minus) sign corresponds to the longitudinal (transverse) coefficient. The coefficients 
are symmetric under particle exchange 1^2, as they should be. Comparing Eqs. [TH] and (THl we 
see that the finite-size effect sets in for r < (a^K^^)^/"^; hence the domains of validity stated in 
Eqs.IUHlll 

Figure [1] shows the coupling diffusion coefficients as a function of inter-particle distance, along 
with their asymptotes. 




Figure 1: Longitudinal (A) and transverse (B) coupling diffusion coefficients as a function of inter-protein 
distance. The diffusion coefficients are scaled by / (Arj^) and the distance by the Saffman-Delbriick length 
K~^. The results for r ^ (a^K~^y^^ (Eq.[T6l solid lines) are presented along with their asymptotes for r ^ 
(Eq. [iTl dotted lines) and {a^K^^y/^ <^ r <^ (Eq. [181 dashed lines). The insets focus on the region of 
shorter distances, where corrections due to protein size become significant (Eq. dash-dotted lines, taking 
m = a2 = 10~^K~^). 



Effective viscosity 

We would like now to calculate the change in membrane viscosity due to the presence of many 
embedded proteins, to leading (linear) order in protein concentration. The definition of effective 



Correlated Diffusion of Membrane Proteins 



7 



viscosity is a subtle point, to which we return in the Discussion. Here we define it using the large- 
distance flow response of the membrane, i.e., we extract the effective viscosity from G(r — )■ oo) —>■ 
G"^^, as it is modified by protein concentration. A similar procedure for a 3D suspension of hard 



spheres (|26f ) correctly reproduces the known effective viscosity of that system, Eq. [5] 



We begin again by applying a point force f at the origin. The membrane flow velocity at position 
r is then given by vf^\r) = Gij{r)fj, where G is given by Eq. [S] Let a disk-like protein of radius a 
be positioned at r'. Due to its finite size it will perturb the membrane flow. The particle is force- 
and torque-free, and, hence, the leading moment of that perturbation is a force dipole, S(r'). From 
our second Faxen-like law, Eq. [T31 we have 

S,j{r') = 2TTfj^a^[d,Gjk{r') + djG,k{r')]fk, (20) 

where we have neglected terms of order a^. This force dipole, located at r', changes the flow velocity 
at position r by Svi{r) = Skj{r')dkGij{r — r'). 

Now suppose that many such mobile proteins are present in the membrane, occupying a fraction 
(j) of its area. The theory being linear, we can superimpose their individual perturbations and 
average over all possible positions r'. This yields an average correction to the flow velocity, 

{Sv,{r))= j dv'p{v')Sk,{v')dkG,j{v-v'), (21) 

where p(r') is the probability density of finding a particle at r'. To leading order in (p we may 
assume a uniform probability density, p{v') = 4)/{'Ka?'). The convolution in Eq. [5T]is then conve- 
niently handled in Fourier space, {5vi[c{)) = [<j) / [TTa'^)\Skj{(\)iqkGij{q) . Substituting Eq.[7]and the 
transform of Eq. [20l we find 

^,(q)=i)f) + (<5i),)=G?ff(q)/,, G<=«=(i-^)g. (22) 

In the limit k — > Eq. [22] reduces to the effective response of a 2D suspension of hard disks, Eq. [6] 
('Sil). For finite k, because of the g-dependent prefactor in Eq. [211 it does not seem at first as if the 
mobile particles could lead to such a straightforward renormalization of the membrane response, as 
they do in 2D and 3D suspensions (|26h. Only in the limit g ^ cxo do we simply get G''*^ (1 — 20)G 
as in the 2D case. In the opposite limit, g ^ 0, we obtain G°* — > G, i.e., the membrane response 
becomes unaffected by the presence of proteins. A closer inspection reveals, however, that Eq. [22| 
could be also obtained from G of Eq. [7| by the following simple substitution (up to linear order in 

</>): 

= G|,,^^^,j^^(i_|.20),k^k(i-20)- (23) 
Furthermore, since k = 2?7t/?7m7 the two substitutions (to linear order in (p) are one and the same. 
Thus, as in a 2D suspension of hard disks (|28|), one can write the effective viscosity of the membrane 
as 

TySf = 77,^(1 + 2<t>), (24) 

provided that this modification is applied to the Saffman-Delbriick length as well. 

One should not be confused by the similarity of Eqs. [6| and [24l the effective hydrodynamic 
response of a protein-laden membrane is not at all similar to that of a 2D suspension. At short 
distances (yet still much larger than the protein size), substituting 77^ — > 7/m(l + 20) and n — > 
k(1 - 2(1)) in Eq. [9] leads to 

r « K-^ : G?7(r) (1 - 2</.)G,, (r) + (25) 

Thus, even in this region of 2D-like behavior, there is an extra term in the membrane response, 
arising from the modification of k. The effect becomes much more dramatic in the large-distance 
limit, where we have from Eq. [10] 

r»«;-i : G?f(r)^G.,(r), (26) 
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without any effect of the embedded proteins. As in Eq. [TOlthe underlying physics is that over large 
distances stresses are transmitted through the surrounding liquid. Hence, the response becomes 
indifferent to the properties of the membrane, be it with or without proteins. 



Concentration corrections to pair diffusion 

The results of the preceding section can be readily used to obtain the corrections to the coupling 
diffusion coefficients due to the presence of many mobile, disk-like proteins, occupying an area 
fraction (j). All we need to do is substitute in Eq.[T6]?7in — > '?m(l + 20), n — > k(1 — 2(/)), and expand 
to linear order in ch. This results in 



SDlir) 



277m K7 



Kr 2 irynry 



1 



'KKr{Kr + 1) 



(27) 



Equation [27| gives the concentration corrections to the bare coupling diffusion coefficients given 
in Eq. [TH Their spatial dependencies are depicted in Fig. [2l The correction to the longitudinal 
coupling is always negative, whereas the correction to the transverse one becomes positive for 
r > . This is because at such large distances the bare transverse coefficient D'^ (Eq. [T7|) 
increases, rather than decreases, with rjm- 




4 

Kr 



-0.1 



B 



I 

4 

Kr 



Figure 2: Concentration corrections to the longitudinal (A) and transverse (B) coupling diffusion coefficients 
as a function of inter-protein distance. The corrections to the diffusion coefficients are scaled by k-^T (j) / {2rji-n) 
and the distance by the Saffman-Delbriick length The results for r 3> (o^k"^)^/"^ (Eq. [27l solid lines) are 
presented along with their asymptotes for r » (Eq. [251 dotted Hues) and (a^K~^)^/'^ <C r- <C k,~^ (Eq. 
dashed lines). 



At large distances we get the following corrections to Eq. [iTl 



r ^ K 



-1 . 



SD- 



7r77ni(Kr)2 



(28) 



Recall that the bare longitudinal coefficient decays at large distances as 1/r, whereas the transverse 
one decays as (Eq.[T7|). Hence, according to Eq.[28l the asymptotic behavior of Z3£ is unaffected 
by the presence of the proteins, D?. This is because at such distances the bare D? is 
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independent of membrane viscosity (Eq. I17p . By contrast, the transverse coefficient is affected 
(increased) by 4> at large distances, — ^ (1 + 2(I))D^. 
At distances much shorter than k^^ Eq. [27] becomes 

(a2K-i)i/3^r<K-i : ,5D£ t W - '/'^^ [ln(W2) + 7 + 1 T 1/2] • (29) 

27r?7m 

Comparing with Eq. I18i we find that the leading term [~ ln(/cr)] is renormalized by (j) as in a 2D 
liquid, D£ rp — > (1 — 20)D£ rp. Yet, the next-order term [0(1)] does not follow the same law. This 
is again because there is another dependence on rjm in the length k^^. 



Discussion 

The theory presented here has yielded several predictions concerning the correlated Brownian mo- 
tion of pairs of membrane proteins, which can be directly checked in two-particle tracking experi- 
ments using Eq. 1151 Equation llGI gives the longitudinal and transverse coupling diffusion coefficients 
as a function of inter-protein distance. An equivalent result, in the form of in-planc response func- 
tions, has been previously reported (jlTl). It is valid at sufficiently large distances and insensitive 
to the size and shape of the proteins. For smaller distances, yet still larger than the protein size 
a, we have derived expressions for the coupling diffusion coefficients to leading order in a/r and 
assuming a disk- like shape of the embedded proteins (Eq. [T9|) . Since the finite-size correction is 
only quadratic in a/r, it becomes significant only for distances r < (o^k^^Y^^. For a ~ 1 nm and 
~ 10'^ nm this crossover length is only ~ 10 nm. Thus, we expect the finite size of proteins 
to affect their hydrodynamic coupling only at nanometer-scale distances. We have provided the 
corrections to the two coupling coefficients due to the presence of other membrane proteins, to 
leading order in their area fraction cj) (Eq. I27p . 

Several particular predictions are worth emphasizing. First, the longitudinal diffusion coefficient 
decays with distance more slowly than the transverse one. Asymptotically, Df^ decays as 1/r, 
whereas decays as fEq. [T7|) . Second, at such large distances Df^ becomes independent 

of membrane viscosity and is, in fact, identical to the longitudinal coefficient in an unbounded 
liquid. This is because the coupling in this regime is mediated by flows in the surrounding liquid. 
Thus, the large-distance longitudinal coupling should be the same for different membranes in the 
same solvent and can be tuned by changing the solvent viscosity. The dominance of stresses 
in the outer liquid holds for the transverse coupling as well, yet this coupling arises from an 
effective force dipole, which is proportional to k"^ = rj^/{2rj[) and, therefore, remains membrane- 
dependent. We note that these asymptotes hold for r ^ k^^ ^ 0.1-1 fim and, hence, may be hard 
to observe in practice. Nonetheless, the difference in the spatial decays should be seen already 
at much shorter distances. (See Figs. [T]4. and [TJ3.) Third, the longitudinal coefficient at large 
distances is predicted to be independent of protein concentration (Eq. [22]). This is merely another 
consequence of the membrane-independence of this coefficient. By contrast, the large-distance 
transverse coefficient not only depends on protein concentration but increases with (Eq. [28] and 
Fig. [213). This unusual result — hydrodynamic interaction enhanced by particle concentration — 
stems from the aforementioned effective force dipole, which increases with ~ 77,11. 

The concept of effective viscosity may have different, not necessarily equivalent, definitions (jlll ). 
The case of a membrane with mobile inclusions seems to clearly demonstrate this difficulty. If one 
were to measure rj^ from the large- distance longitudinal coupling coefficient as a function of (p, 
one would find no concentration effect, yyjjf = rj^. If, however, rj^ were extracted from the large- 
distance transverse coefficient, a concentration dependence identical to that in a 2D liquid would 
follow, 77^^ = (1 -(- 2(j))r]^. (Compare Eqs. [T7]and[25]) The calculation by Henle and Levine (jsoh 
is based on the same (Saffman) theory, yet follows Einstein's original definition of the effective 
viscosity as the coefRcient relating average stress (or dissipation rate) with strain rate under a 
given global shear ffow |27|). They find yet another concentration dependence in the limit Ka — > 0, 



Correlated Diffusion of Membrane Proteins 



Vm = ??m(l + 30). In the cases of 3D and 2D suspensions, as shown, respectively, in Ref. (|26ll and 
the current work, the aforementioned three definitions of yy^f give identical results. The different 
behavior of membranes lies in the fact that they do not conserve momentum, i.e., in the appearance 
of the length scale and its dependence on rj^^. 

Thus, the effect of inclusions on membrane viscosity depends on the definition of that transport 
coefficient and the experiment under consideration. The definition used here relates to how the 
response of the membrane to local perturbations changes with 4>. It is relevant, therefore, to 
microrheological and particle-tracking experiments. Einstein's definition, as used by Henle and 
Levine (|30[). should be appropriate for larger-scale rheological measurements. We have found that 
the effective response in the presence of proteins is given, as in a 2D suspension, by substituting 
Vm = '7m(l + 2(/)). Yet, unlike 2D suspensions, this substitution should be made also in 

the Saffman-Delbriick length, k^^ = rin-i/{2rii) ?7°f /(2?7f). The extra dependence of k on leads 
to a qualitatively different effective response. A clear demonstration is given by the concentration 
correction to the transverse coupling, SD^, where the interplay between the (/)-dependencies of Tym 
and K leads to a sign reversal of that term (Fig. 03). Another consequence is the aforementioned 
independence of the large- distance longitudinal coupling, Z?£, on protein concentration. 

We note that a similar absence of a concentration effect on the large-distance response has been 



observed in another quasi-2D system — a suspension confined between two plates (|3ll |32| ). The 



physical origins of the two phenomena, however, are slightly different. In the confined suspension 
momentum is lost to the solid boundaries, and the far response arises solely from liquid mass 
displacement, which is not affected by the presence of particles. In membranes the far response does 
arise from momentum diffusion, yet these dynamics take place in the outer liquid and, therefore, 
are insensitive to the presence of membrane inclusions. 

The current analysis has involved several rather severe approximations. It should be regarded, 
therefore, as a first step toward understanding the correlated dynamics of membrane proteins or, 
alternatively, as a possible means to isolate simple hydrodynamic effects from other, more specific 
ones. First, we have focused on the hydrodynamic coupling between proteins, neglecting any 
direct interaction (|33l). Such interactions may arise from actual (e.g., electrostatic) potentials or be 
induced by the perturbations that inclusions introduce in the membrane (fiol. Is^. lisl). Second, as 
in previous theories, we have considered a homogeneous membrane, whereas actual biomembranes 
are believed to contain various heterogeneities and domains (|3)- A homogeneous hydrodynamic 
description may still be applicable inside such a sub-micron domain. Third, the calculations have 
been made in the limit of very small particle size, Ka <^ 1. As the typical values of are micron- 
scale, this should be a good approximation for practically all membrane proteins. The Saffman 
theory can be extended to large values of Ka as well, yet the calculations become significantly more 
complicated (1^ 13) ■ Fourth, we have treated the membrane as a perfectly flat surface, whereas 
in practice it is curved and fluctuating. Curvature and bending fluctuations, apart from their 
aforementioned ability to induce interactions between embedded proteins, may also affect their 
2D-projected diffusion as observed in experiments (37, SSh . 

Finally, we have studied the effect of protein concentration to linear order only. As in the 
much simpler case of a 3D suspension of hard spheres, extension to higher orders in (/) should 
be difficult, involving static and dynamic correlations between particles. Nevertheless, some of 
our results clearly emanate from more fundamental considerations and, therefore, should hold for 
higher values of (p as well. For example, the (/<- independence of I?£(r ^ k^^) arises from stresses 
being transmitted through the outer liquid; thus, we conjecture that it is valid to all orders in (f>. 
It is also plausible that, in the limit /ca <C 1, expressions of higher order in (p could be obtained by 
merely substituting for rym (both directly and in k) the effective viscosity, r;^ (</)), as calculated for 
a 2D suspension. 
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Appendix A: Faxen laws for a membrane 

In this appendix we calculate the approximate membrane-analogue of the 3D Faxen laws. These 
laws relate the linear velocity U and angular velocity H of a rigid particle, and the flow field v(r) 
in which it is embedded, with the moments of force distribution that it exerts on the embedding 
fluid. 

Let the center of a disk of radius a be located at the origin, and let f(r') be the distribution 
of forces that it exerts on the system. This force distribution changes the membrane velocity at 
position r by /c?r'Gy (r — r')/j(r'), where G is the membrane-analogue of the Oseen tensor, Eq. 
[8l Assuming no slip at the particle perimeter, we have 



U, + {flx r), = i-,(r) + J dr'Gy(r - r')/j(r'). 



(30) 



We have intentionally left the domain r' of the force distribution unspecified, since the calculation 
is insensitive to it; the vector r, however, must be on the particle perimeter, where the no-slip 
boundary condition is imposed. 

Assuming that v(r) changes very moderately on the scale of a, we expand v(r) ~ v(0) + 
afidi\'{0) + ^a^fifjdidjw{0). Since both r and r' in Eq. [30] are of order a, and we have been 
assuming throughout this work <^ 1, we can substitute for G its short-distance asymptote, 
Eq. [51 Within these approximations, integrating Eq. [301 over r yields the analogue of the first 
Faxen law, 



47r?7n 



ln(Ka/2) + 7 + 0(Ka)2 



1 



v(0) + -a^V^v(O) - U + 0{a^V^^^ 



(31) 



where F = J dr'f(r'). 

Next we multiply both sides of Eq.[3D|by r and integrate over r. Separating the resulting tensors 
into symmetric and antisymmetric contributions, we obtain the analogue of the second Faxen law 
for the torque and force dipole (stresslet). 



Sij — 2nri^a [1 + 0{Ka) 



l + ia2v2 + 0(a4v4) 
8 



[V X v(0)] - 2nj (32) 
[d^v,{0) + dMO)], (33) 



where L = / dr'[r' x f (r')] and S^j = (1/2) / dr'[r'Jj{r') + r;./,(r')]. 

In regular Stokes flows k = 0, V^v — 0, V^(V x v) = 0, and the Faxen laws become exact. Thus, 
Egs. [3lH331 are exact for 2D liquids, where, additionally, the term proportional to V^(V x v) in 
Eq. 1321 vanishes. For membranes, however, their validity is restricted to sufficiently small particles, 
Kfl <^ 1. In addition, the terms of order are valid provided that the considered flow is not too 
uniform, \V'^v/v\ 3> k^. 
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